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Computer  Simulations  of  Microwave  Generation 
From  Relativistic  Electron  Beams 


This  is  the  Final  Technical  Report  on  work  performed  under  the  support  of  the  Air 
Force  Office  of  Scientific  Research  on  the  subject  of  simulating  microwave  generation. 
This  report  contains  the  write-ups  that  we  sent  for  the  upcoming  book  on  power  vacuum 
electronics.  The  first  section  (Chapter  10)  describes  the  simulation  algorithm  based  on 
mode  expansion  technique  developed  at  UCLA  while  the  second  section  (Chapter  5) 
gives  the  real  amplification  mechanism  for  the  Phigtron. 


I.  Mode  Expansion  Technique  (METPIC) 

A.  Introduction 

Using  the  Coulomb  gauge,  an  algorithm  for  solving  Maxwell’s  equations  that  is 
capable  of  isolating  the  solenoidal  current  component  from  the  total  current  has  been 
developed.  The  electromagnetic  field  arising  from  the  solenoidal  current  is  expanded  in 
terms  of  Bessel-Fourier’s  functions.  It  is  sufficient  to  retain  only  a  few  dominant  global 
electromagnetic  modes  to  simulate  most  microwave  devices.  Electrostatic  modes  that  are 
localized  within  the  beam  region  are  followed  by  utilizing  B-splines.  This  mode 
expansion  technique  is  most  effective  in  addressing  basic  physical  issues  involved  in 
microwave  devices  with  simple  boundary  conditions. 

The  electron-wave  interaction  in  microwave  devices  is  a  highly  nonlinear 
phenomenon  and  requires  computer  simulations  to  follow  its  evolution.  The  electron 
convection  current  calculated  in  simulation  often  contains  both  the  electromagnetic  and 
electrostatic  components.  This  total  current  is  usually  used  as  the  source  for  exciting  the 
electric  and  magnetic  fields  in  electrodynamic  structures.  With  this  approach,  it  is 
important  to  make  sure  that  the  equation  of  continuity  is  satisfied  at  all  times.  In 
MAGIC1,  this  is  accomplished  by  correcting  the  electric  field  at  every  time  step  to  ensure 
that  Poisson's  equation  is  satisfied.  Most  of  the  other  approaches  employed  to  simulate 
microwave  devices  are  for  medium  electron  beam  current  and  assume  fixed  transverse 
field  profiles  determined  from  the  boundary  conditions  of  the  structure.  The  electrostatic 
field  evaluated  from  the  Poisson  equation  is  often  added  to  the  fields  generated  by  the 
total  current.  In  doing  so,  the  dominant  mode  of  space  charge  fields  may  be  doubly 
counted.  In  order  to  be  able  to  correctly  investigate  the  effects  of  space  charge  field  on 
the  performance  of  microwaves  devices,  a  new  algorithm  for  solving  Maxwell's 
equations  which  is  capable  of  isolating  the  solenoidal  current  from  the  total  current  have 
been  developed. 

Since  early  seventies,  particle  simulations  based  on  the  Fast  Fourier  Transform 
(FFT)  technique  to  solve  Maxwell’s  equations  have  been  in  use  in  plasma  physics  and 
radiation  generation2,3  at  Princeton  Plasma  Lab  and  UCLA.  In  1982,  Caplan  applied  the 
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mode  expansion  technique  to  simulate  Gyrotron  Devices4.  In  microwave  devices,  a 
cylindrical  waveguide  is  often  provided  to  confine  the  electromagnetic  wave  propagation 
along  the  axial  direction.  A  combination  of  Bessel's  functions  can  always  be  found  to 
satisfy  the  proper  radial  boundary  conditions  and  Maxwell's  equations.  The  field  quantity 
can  be  expanded  into  a  series  of  Bessel's  functions  that  are  complete  orthogonal 
functions.  For  electromagnetic  wave  problems,  the  field  in  the  exterior  of  the  electron 
beam  takes  the  form  of  a  wave  moving  outwards  if  there  is  no  structure  to  confine  the 
wave  propagation.  In  a  cylindrical  waveguide,  the  transverse  mode  becomes  a  global 
mode  and  retaining  a  few  modes  is  often  sufficient  to  evaluate  the  performance  of 
microwave  devices.  Along  the  axial  direction,  the  Fourier  series  expansion  can  be  used 
to  represent  the  self-consistent  field  profile  that  plays  an  important  role  in  determining 
the  efficiency  of  microwave  devices. 

B.  Solenoidal  Current  Component 

By  using  the  Coulomb  gauge5,  the  current  density  (or  electric  field)  can  be 
decomposed  into  two  terms, 

t-lt+if  .  0) 

Ji  is  the  longitudinal  (irrotational,  V  x  it  =0)  current  component  and  is  the  driving 

source  for  electrostatic  fields,  while  Jt  is  the  transverse  (solenoidal,  V  •  Jt  =0)  current 
component  and  is  the  driving  source  for  electromagnetic  fields.  The  magnetic  field,  of 
course,  only  has  transverse  component.  Starting  from  the  vector  identity, 


Vx(Vxf)  =  V(V*t’)-V2lf  , 
the  transverse  current  can  be  determined  from  the  total  current 


(2) 


V2(jf)  =  -VxVxT  =V2?  -V(V»?  )  . 


The  solution  of  Eq.  (3)  in  integral  form  is5 


it  =~  VxVx 
47t 
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("x  -?) 


dV  . 
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This  form  of  solution  is  difficult  to  obtain  numerically.  In  the  differential 
formulation,  Eq.  (3)  can  be  easily  solved  if  the  field  quantity  is  expanded  into  a  series  of 
complete  orthogonal  functions.  For  example,  sets  of  such  functions  are  the  Fourier  series 
and  the  Bessel  functions  which  are  the  solutions  for  cylindrical  waveguides. 

Consider  an  infinite  and  periodic  system.  In  this  case,  it  is  convenient  to  solve 
Maxwell's  equations  in  the  Fourier  space3  and  the  splitting  of  the  current  into  transverse 
and  longitudinal  components  can  be  easily  done.  From  Eq.  (3),  the  transverse  current  can 
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be  written  as 


(5) 

where  is  the  wave  vector.  Based  on  this  approach,  for  the  last  three  decades, 
numerous  interesting  nonlinear  phenomena  in  plasma  physics3  have  been  unraveled. 

In  microwave  devices,  a  cylindrical  waveguide  (r,0,z)  is  often  employed  to  guide 
the  electromagnetic  wave  axial  propagation.  The  electromagnetic  fields  in  a  cylindrical 
waveguide  with  a  constant  cross-section  can  be  expanded  in  terms  of  TE  and  TM  modes 
as  follows.  For  TE  modes4 


t_=  2  E^fct)  ( £  x  v.cl) 

p=l,2 


mn 
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^mn  =  -  Z  [BP  (z,t)V,CP  +  z  B  P  kmnCP  ], 

mn  L  Tmnv  ’  '  1  ran  Lmn  mn  J 


P=l,2 


where  the  circular  wave  functions  for  the  two  polarizations  are  given  by 


C1,2  =  Jm  [  ]  [  cos  m9  ,  sin  m0  ] 

mn  Kmn  i\v 
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and  J'm  (k  mn)  =  0.  Here,  rw  is  the  waveguide  radius.  For  TM  modes 


^mn  -  -  Z  [E-rn^2’1)  VlCnm  +  Z  ELPmn  kmn  C™  ]  ’ 

p=l,2 


mn 


mn 
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B  mn  Z  B  (Z,t)  (  Z  X  VlCJ  , 

P=l,2 

where  Jm(  kmn )  =  0. 

Take  the  spatial  variation  of  the  electron  beam  current  to  conform  the  electric 

field  variation.  For  TE  modes,  substituting  Eq.  (6)  into  Eq.  (3)  yields  V(V*~?)  =  0 

and  thus  =~t  (excluding  the  spatially  uniform  current) .  For  TM  modes,  as  a  simple 
example,  an  algorithm  under  the  assumption  of  axial  symmetry  is  developed.  In  this 
situation,  the  electric  fields  are  along  the  z  and  r  directions  and  only  these  current 
components  (Jz  and  Jr)  needs  to  be  retained  as  the  source  for  exciting  the  electromagnetic 
field 


where  kon  is  the  zero  of  the  Bessel  function  Jo. 

The  second  term  of  the  right-hand  side  of  Eq.  (3)  becomes 


rA,  .  2  T  ,  ikzkon  T  '  .  .K0nr 

V(V«J)  =  [z(-kzJz  +  Jr)  Jo(  r  ) 

1 W  UV 
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The  solenoidal  current  (Eq.  3)  can  thus  be  written  as 


2  2 
r  t  kp»  T  ,  ikzkpn  T  \  T  /  i_  _L  \  _ikzz  //  kpn  ,  2i 
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Notice  that  Eq.(l  1)  is  divergence  free.  These  currents  are  retained  as  the  source  to 
advance  the  electromagnetic  fields  through  Maxwell's  equations. 


C.  Nonuniform  Waveguide  Wall 

The  nonuniformity  in  waveguide  wall  causes  the  modes  of  the  waveguide  to 
couple  among  each  other.  There  are  at  least  two  different  ways  to  describe  this  mode 
coupling  process.  In  the  configuration  that  uses  a  taper  to  connect  two  uniform 
waveguide  such  as  the  extended  interaction  cavity,  the  most  convenient  way  is  to  expand 
the  electromagnetic  field  of  the  actual  waveguide  in  terms  of  so  called  local  normal 
modes.  These  modes  are  identical  in  form  with  the  exact  mode  solutions  of  the  uniform 
guide.  However,  the  guide  dimensions  appearing  as  parameters  in  the  mode  expressions 
and  the  eigenvalue  equation  are  now  allowed  to  be  functions  of  z  to  conform  exactly  to 
the  actual  shape  of  the  nonuniform  waveguide.  These  modes  are  not  themselves 
solutions  of  Maxwell’s  equations  because  their  parameters  are  function  of  z.  However, 
they  can  be  superimposed  to  yield  the  desired  solution.  This  approach  results  in  mode 
coupling  coefficients  that  vanish  everywhere  except  on  the  tapered  portion  and  is  used  in 
simulating  the  extended  interaction  cavity.  The  other  way  is  to  use  the  normal  modes  in  a 
uniform  waveguide  as  the  base  for  field  expansion.  This  latter  approach  is  more 
convenient  for  handling  waveguide  with  periodic  varying  cross-section. 

Substituting  Eqs.  (6)  and  (7)  into  Maxwell's  equations  in  three  dimensions,  for  each 
TE  mode,  their  governing  equations  can  be  reduced  to  an  equivalent  set  of  one 
dimensional,  time  dependent,  inhomogeneous  transmission  line  equations4,6. 
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Here  cocmn  =  ckmn/rw(z)  is  the  waveguide  cutoff  frequency, 

( —  f  =  5'87J  X2~~ ,  K±02=  rw(°c-—  ,  Umn  =  Jm2(kmn)  (  1  -  m2/kmn2 ) 

03C0  K_j_o  pz0  C 

rr  _  ^  kmn  Jm  (  kml  ) 

mnl "  (kml2-kmn2)Jm(kmn) 

and  the  electron  beam  is  .modeled  by  a  collection  of  Np  macroelectrons  distributed  over 
length  Lp.  Each  macroelectron  represents  N*Lp/Np  electrons  where  N  is  the  number  of 

electrons  per  unit  length.  Here,  c  is  the  speed  of  light,  coco  is  ckmn/rw(0),  Vj  is  the 
electron  velocity,  and  pz0  is  the  average  beam  axial  velocity  divided  by  c.  In  order  to 
avoid  the  singularity  at  cavity  axis  in  advancing  electron  trajectories,  the  calculations  are 
carried  out  using  Cartesian  coordinates.  The  x  and  y  components  of  the  transverse  field 

functions  exmn  and  bimn  are  given  by 
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The  electromagnetic  fields  for  the  TEmn  mode  acting  on  the  electrons  are  then  given  by 
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Following  the  approach  of  Solymar7,  the  waveguide  radius  in  this  scheme  is  allowed  to 
be  an  arbitrary  continuous  function  of  axial  position,  allowing  one  to  model  profiled 
interaction  circuits  and  mode  conversion.  The  axial  spatial  dependence  of  the 
electromagnetic  field  is  still  represented  by  the  Fourier  series  expansion  and  the  leap-frog 
scheme  for  electron  pushing  in  PIC  codes  is  retained.  The  out-going  wave  boundary 
conditions  are  imposed  along  the  axial  direction.  The  code  is  stretched  in  the  sense  that 
the  electromagnetic  field  allows  only  one  dimensional  (axial)  spatial  evolution  but  is  used 
to  simulate  three  spatial  dimensional  phenomena.  This  version  of  code  is  called  stretched 
PIC  code  in  plasma  physics.  These  sets  of  equations  appeared  first  in  1982  in  Caplan’s 
Ph.D  dissertion3.  He  also  used  an  effective  damping  rate  for  the  electric  component  of 
the  electromagnetic  field  to  simulate  the  waveguide  wall  resistivity. 


D.  Space  Charge  Field 

In  a  cylindrical  waveguide,  the  transverse  mode  is  a  global  mode  and  a  finite 
number  of  modes  is  often  sufficient  to  represent  the  electron-electromagnetic  wave 
interaction.  On  the  other  hand,  for  electrostatic  problems  the  fields  outside  the  beam 
location  decay  exponentially,  with  no  spatial  oscillation  and  no  wave  propagation.  To 
simulate  this  localized  field,  a  great  many  transverse  modes  are  needed  along  the  radial 
direction.  As  a  result,  the  scheme  of  using  finite-differencing  or  B-spline  are  a  better 
approach  to  evaluate  the  electrostatic  field.  The  B-spline8'10  of  fourth  order  will  be  used 
to  solve  the  two-dimensional  Poisson  equation  in  cylindrical  coordinates 


a2^  1  so  a2® 

~S?  +7i;  +_a?  =-4nP  • 


(15) 


with  the  boundary  conditions  of  O(r,0)  =  <T(r,L)  =  <I>(rw,z)  =  0. 
and  p  along  the  z-direction  in  terms  of  sine  functions 
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Expand  the  functions  <b 


z  ,  (16) 


where  Nz  is  the  number  of  grid  points  in  the  z-direction.  Both  sine  and  cosine  functions 
should  be  used  in  the  case  of  periodic  axial  boundary  conditions.  Substituting  Eq.(16) 
into  Eq.(15)  yields 


1  k7T 

#k"  +  7  Ok'-(-j-)2  0>k  =  - 47tpk 


(17) 


where  prime  denotes  taking  derivative  with  respect  to  the  r  coordinate. 

The  collocation  method  with  B-splines  is  used  to  numerically  solve  Eq.(17).  Let 

{  n  }^r  be  the  radial  grid  points  (knots),  0  =  r,  <  ....<  rNf  =  rw.  The  solution  of  Eq.(17) 
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can  be  approximated  by  a  spline  S  of  fourth  order  on  the  grid  {rj } .  A  spline  can  be 
expressed  in  the  following  form 

n 

S(r)  =  ZaiBi(r)  ,  (18) 

i=l 

where  Bj(r)  are  so  called  B-splines  which  are  basic  functions  of  the  space  of  all  cubic 
splines  on  the  grid  {p} .  n  =  Nr  +  2  is  the  dimension  of  this  space. 

A  Spline  function  is  a  function  consisting  of  polynomial  pieces  on  subintervals, 
joined  together  with  certain  smoothness  conditions  and  is  used  in  PIC  code  for 
interpolation. 


o  (  1  if  ti  <  r  <  ti+1 
Bi  (r)  =  { 

0  otherwise 

This  is  the  B-splines  of  degree  0.  To  generate  all  the  higher  B-splines,  the  simple 
recursive  definition  can  be  used 

k  r-ti  k_l  ti+k+i-r 

B-w-W*  <r)  +  (wft^)Bi*'(r)  tel) 

n 

Bi^(r)  =  1  for  all  r 

i=l 


Substituting  S(r)  for  O(r)  in  Eq.(16)  yields 


V  (k)rT.»  1  '  .  k7T  2  n  4 

2-i  ai  [Bi  +-  Bi  -(  —  )  Bj  ]  =  - 4Ttpk 


i=l 


(19) 


Since  the  values  of  pk  are  known  at  the  grid  point  rj,  we  will  evaluate  Eq.(19)  at  these 
points  to  determine  the  coefficients  aj(k). 


E(k)  it  1  f  krc  ^ 

ai  [  Bj  (rj)  +  -  Bj  (q)  -  ( — )  Bi(rj)  ]  =  -  4Tipk(rj)  , 
i=l  J 

j  =  2,  3, ...,  Nr  (20) 

O' 

Because  the  term  —  in  Eq.(ll)  is  singular  at  r=0,  we  must  treat  it  carefully.  The 
property  of  axis-symmetry  implies  Er=-O-0  at  r=0.  The  L'Hospital  rule  is  used  to 
evalute  this  singular  term 
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(21) 


lim  -  O'  =  lim  O"  =  0"(  0  )  . 
r— >0  r  r->0 

Using  Eq.(21),  Eq.(20)  at  r=0  becomes 

n 

Z  ai(k)  [  2Bj"(0)  -  ( Y )  Bj(0)  ]  =  -  4rcpk(0).  (22) 

i=l 

Besides  Eq.(20)  and  (22),  we  have  two  additional  boundary  conditions  O'(O)  =  0  and 
Ok(rw)=  0.  In  terms  of  B-splines,  these  become 
n  n 

Z(k)  (k) 

aj  Bj'(0)  =  0  and  2,  aj  Bf(rw)  =  0. 

i=l  i=l 

00 

Now  we  have  altogether  Nr+2  equations  for  Nr+2  unknown  coefficients  a,  and  the 
resulting  equation  can  be  written  in  the  form  of  tridiagonal  matrix.  After  aj(k)  have  been 

evaluated,  the  electrostatic  field  can  be  determined  from  E^  =  -VO.  In  comparison  with 
the  conventional  finite  difference  scheme,  the  present  approach  requires  less  storage 
spaces  and  is  more  accurate. 


E.  Simulation  Results 

In  order  to  validate  the  PIC  code  based  on  the  mode  expansion  techniques,  simulation 
results  were  compared  with  experiment,  MAGIC,  and  theoretical  prediction.  In  the 
following,  only  crucial  results  for  the  purpose  of  comparison  are  presented,  the  other 
detailed  can  be  found  in  their  respective  references. 

a.  Nonlinear  Mode  Interaction 

In  the  experimental  result  on  gyrotron  traveling  wave  amplifier  operated  at  the 
TEn  mode,  it  is  observed  that  by  increasing  the  injected  power  of  the  TEn  mode,  the 
TE2i  mode  self-oscillation  can  be  suppressed(Fig.la).  The  following  parameters  were 
used  in  the  experiment  and  the  simulation:  Vb=87.75kV,  Ib=T.25A,  a=l,  rw=0. 2654cm, 
B0=12.38kG,  and  f0=34.7  GHz.  Since  low  current  beam  was  used,  only  the  TE  modes 
were  retained  in  simulations  (space  charge  field  is  neglected)  and  the  result  was  shown  in 
Fig.  lb.  The  agreement  between  them  are  very  impressive.  As  a  result,  mutual 
verifications  between  the  simulation  and  experiment  were  established  and  have 
strengthened  the  understanding  of  the  multimode  behavior  and  added  to  the  predictability 
of  high-power  gyro-TWT's.  Much  physical  insight  into  the  multimode  dynamics  has 
been  obtained  with  the  valuable  tool  of  numerical  diagnostics  which  can  be  easily 
performed  in  simulations.  The  mechanisms  for  mode  competition  /suppression,  for 
example,  are  clearly  seen  through  the  diagnostics  of  the  perturbed  beam  distribution  in 
momentum  and  energy  spaces.  Simulation  results  (Ref.[ll])  clearly  show  that  the 
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organized  beam  perturbation  associated 
3  with  one  mode  appears  as  energy  and 
velocity  spread  to  another  mode  of 

2  J5  different  frequency  or  mode  structure,  and 
^  vice  versa. 

1  Mode  competition  processes  are 
examined  from  a  different  perspective  in 

o  Fig.  lc,  which  plots  the  temporal  behavior 

O  S 00  1OQ0  1500  2000  o  ?  r  i 

of  the  spatially  integrated  wave  power(P0) 
of  the  TEn  and  TE21  modes.  As  soon  as 

3  the  electron  beam  is  turned  on  at  t  =  0,  the 
TE21  mode  grows  spontaneously  to  a 

2  jg  steady  state  following  a  few  bounces 
«  caused  by  the  instant  turn-on  of  the  beam. 

1  At  t  =  42.5  nsec,  a  900  W  drive  is  injected, 
which  grows  to  saturation  in  ~  2  nsec,  a 
o  time  scale  comparable  to  the  beam  transit 

O  500  IOOO  1500  2000  .  .  ,  -  .  . 

time  through  the  waveguide.  This  is 
expected  since  the  beam  velocity  and  the 
TEn  wave  group  velocity  are 
approximately  equal  under  the  grazing 
resonance  condition.  Suppression  of  the 
TE21  oscillation  is  mostly  completed  at  t  = 
46  nsec.  But  it  is  not  until  the  TEn  wave 
grows  to  a  significant  level  that  the 
oscillation  level  starts  to  diminish,  a  fact 
reflecting  the  nonlinear  nature  of  mode¬ 
mode  interactions.  As  the  drive  is  turned 
Fig.  1  Nonlinear  mode  interaction  in  off  at  t  =  58  nsec,  TE21  oscillation 

Gyro  amplifier  (a)  Measured  TEn  immediately  comes  back  to  the  pre-drive 

output  power  (filled  circles)  and  TE2i  level.  This  numerical  sequence  has 

oscillation  power  (crosses)  vs.  the  TE]  1  precisely  verified  the  experimental 

drive.  (b)  Simulated  data  of  the  behavior, 
measurements  in  (a),  (c)  Time  evolution 

of  mode  power.  The  beam  is  turned  on  b.  Self-field  Effects  on  High  Power 

at  the  reference  time  t  =  0.  The  900  W  Gyro  Backward  Wave  Oscillators 

TEn  drive  is  turned  on  at  t  =  42.5  nsec 

and  turned  offatt  =  58  nsec.  In  a  recent  MAGIC  simulation12 

carried  out  at  Phillip's  Laboratory,  on  the 
performance  of  high  power  fundamental  gyro  backward  wave  oscillators  with  conducting 
axial  boundary  conditions,  the  following  parameters  were  employed  Vb  =  350  kV,  lb  = 
800  A,  a  =  1.5,  rw  =  4.37  cm,  Bo  =  3  kG,  L  =  15  cm,  and  guiding  center  rg  =  0.44  rw. 

In  this  high  current  device,  the  space  charge  filed  plays  an  important  role  in 
determining  the  efficiency.  To  verify  that  the  present  algorithm  for  evaluating  space 
charge  field  is  accurate,  a  simulation  using  the  same  parameters  that  were  used  in 
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MAGIC  simulation  was  performed.  The 
spatial  axial  momentun  distribution  at 
coct  =  300  from  our  simulations13  that 
included  both  EM  and  ES  fields  (Fig.2a) 
is  compared  with  the  result  from 
MAGIC  code  (Fig.  2b).  The  agreement 
between  the  two  results  using  different 
algorithm  is  quite  well. 


5x10s 


o 


— 5x1 08 
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c.  Plasma-filled  Microwave  Devices 

To  illustrate  the  effectiveness  of 
separating  the  electric  field  into  the 
electromagnetic  and  electrostatic 
components,  a  device  operated  at  the 
TM  mode  must  be  simulated.  In  a  recent 
long  pulse  (x  =  120ps)  PASOTRON 
(Plasma  Assisted  Slow-Wave  Oscillator) 
experiment14,  the  plasma  was  observed 
to  totally  neutralize  the  beam  space 
charge  field  while  no  plasma  return 
current  was  generated.  The  generation 
of  plasma  return  current  can  occur  during  the  beam  tum-on  and  turn-off  periods.  This  is 
because  a  time  varying  beam  current  (Jb)  will  induce  a  time  varying  Be  which  in  turn  will 
produce  an  inductive  Ez  field  to  drive  the  plasma  current.  However,  the  effect  arising 
from  the  edge  of  the  beam  pulse  may  not  be  important  in  a  long  pulse  beam.  The 

remaining  JbZ  x  Be0  force  is  capable  of  pinching  the  electron  beam  and  thus  substitute 
as  the  solenoidal  magnetic  field  to  confine  the  beam  transport.  In  the  experiment,  the 
following  parameters  were  used:  Vb  =  96kV,  lb  =  200A,  average  waveguide  radius  rw  = 
3.25cm,  wall  ripple  amplitude  n  =  ,715cm,  and  ripple  period  zo  =  2.4cm.  The  dispersion 
relation  and  growth  rate  for  a  periodic  system  using  these  parameters  are  given  in  Fig.  3a 


o  a  io  is 


Z<cm) 

Fig.  2  Comparison  of  electron  momentrum 
axial  spatical  distribution  from  (a)  our  code 
and  (b)  MAGIC  code. 


coc 

in  which  the  plasma  density  is  assumed  to  be  <bp  =  0.6coc  (fc=  ^  :  waveguide  cutoff 

frequency).  Due  to  the  difficulty  in  measuring  the  plasma  density  when  the  electron- 
wave  interaction  is  taking  place,  the  information  on  plasma  density  was  not  given  in  the 
experiment.  The  smooth  wall  Trivelpiece-Gould  dispersion  curve  is  a  good 
approximation  for  plasma  waves  if  the  plasma  column  as  a  result  of  pinching  is  far  away 
from  the  rippled  wall  boundary.  This  is  due  to  the  fact  that  plasma  waves  are  localized 
inside  the  plasma  column.  The  intersection  "A"  is  the  slow  electromagnetic  wave  (TMoi) 
and  is  absolutely  unstable.  If  either  beam  current  or  interaction  length  exceeds  the  start 
oscillation  threshold,  the  wave  will  grow  from  initial  noise  fluctuations  and  eventually 
saturate.  On  the  other  hand,  the  intersection  "B"  is  convectively  unstable.  The  wave  will 
grow  from  the  initial  noise  fluctuation  at  the  entrance  and  convect  out  of  the  interaction 
region  without  strongly  interfering  with  the  evolution  of  mode  "A".  However,  the  growth 
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rate  of  mode  "A"  is  rather  small  even  in 
an  infinite  system  (Fig.  3a).  In  a  real 
experiment,  the  finite  length  effect  will 
further  reduce  the  growth  rate.  It  is  quite 
possible  that  without  the  reflection  from 
mismatch  between  device  components  or 
other  effects  caused  by  the  presence  of 
plasmas  to  enhance  the  growth  rate,  the 
oscillator  might  be  stable.  Computer 
simulations  did  demonstrate  that  the 
presence  of  plasma  is  capable  of 
enhancing  the  oscillator  gain  in  addition 
to  neutralizing  the  beam  space  charge 
field.  However,  only  the  result  relevant 
to  the  electric  field  decomposition  is 
presented.  The  frequency  spectrum  of 
the  electromagnetic  wave  determined 
separately  from  beam  and  plasma 

transverse  currents  ( Jt  )  is  shown  in  Fig. 
3b  for  (Op  =  0.6  coc.  The  spectrum 

contains  only  one  peak  at  about  1 .3  coc 
for  both  cases  and  is  very  close  to  the 
frequency  of  intersection  "A". 
Simulation  results  show  no  trace  of 
space  charge  field  and  clearly 
demonstrate  that  the  field  arising  from 
the  solenoidal  current  contains  only  the 
electromagnetic  component.  The 
electrostatic  frequency  spectrum  (Fig. 
3b)  evaluated  from  beam  electrons  peaks 
at  the  same  frequency  as  that  of  the  EM 
wave.  This  is  due  to  the  synchronism 
between  the  beam  slow  space  charge 
wave  and  the  slow  electromagnetic  wave. 


(a) 


u/cjc  6>/cjc 


Fig.  3  Illustration  of  the  electric  field 
separation  into  EM  and  ES  components  using 
a  plasma-filled  device  as  an  example  (a)  The 
dispersion  relation  of  the  TMoi  mode,  beam 
slow  space  charge  wave,  and  plasma 
Trivelpiece-Gould  mode.  The  growth  rate  for 
rb  =  0.25  cm  is  also  shown,  (b)  Frequency 
spectrum  of  EM  and  ES  fields  from  beam  and 
plasma. 

There  is  a  secondary  peak  at  the  frequency  of 


beam  fast  space  charge  wave.  In  the  absence  of  BWO  instability,  the  electrostatic  plasma 
response  peaks  at  the  frequency  of  the  Trivelpiece-Gould  mode.  The  peak  shifts  to  the 
frequency  of  beam  slow  space  charge  wave  when  the  BWO  instability  is  excited.  The 
result  of  frequency  spectra  clearly  demonstrates  that  the  scheme  of  separating  EM  and  ES 


works  nicely. 


d.  Mode  Conversion 


The  inhomogeneity  of  the  waveguide  cross-section  in  the  tapered  section  provides 
the  coupling  (mode  conversion)  among  various  waveguide  modes.  This  mode  coupling 
process  could  affect  the  performance  of  the  extended  interaction  cavity  by  determining 
the  mode  purity  in  electron-wave  interaction  for  each  individual  cavity.  In  order  to  be 


sure  that  the  stretched  PIC  code  takes  into 
account  the  mode  conversion  properly, 
simulation  results15  were  compared  with 
the  theoretical  prediction  of  Ref.  [16]. 
According  to  Eq.  (12),  the  coupling 
strength  is  proportional  to  the  slope  of  the 
taper  and  inversely  proportional  to  (ko/2  - 
kon2).  In  order  to  observe  substantial  mode 
conversion,  a  30°  taper  angle  and 
conversion  between  the  TE02  and  TE01  are 
chosen  as  an  example  for  comparison.  It 
is  straightforward  in  simulation  to  include 
coupling  to  additional  modes  if  necessary. 
The  dimensions  of  gyrotron  oscillator  for 
simulation  are  shown  in  Fig.  4a.  The 
mode  profiles  for  the  TE02  and  TE01  modes 
from  theoretical  analyses  (Ref.  16)  are 
displayed  in  Fig.  4b.  A  70  kV,  10  A  beam 
with  a  =1.5  and  guiding  centers  located  at 
0.096  cm  and  B0=  35.54  kG  were  used  for 
simulations.  The  field  profiles  from  our 
stretched  PIC  code  simulation  are  shown 
in  Fig.  4c  which  agrees  quite  well  with 
Fig.  4b. 


z(mm) 


Fig.  4  Comparison  of  mode  conversion  in 
taper  region,  (a)  Cavity  dimensions,  L2  = 
1.5  cm,  Rw  =  0.37  cm,  0,  =  0.5°,  02  =  30°. 
(b)  Field  profiles  (TE01  and  TE02)  from 
theory,  (c)  Field  profiles  from  simulation. 


II.  Amplification  Mechanism  in  The  Output  Section  Of  The  Harmonic 
Multiplying  Gyrotron  Traveling  Wave  Amplifier 

Particle-in-cell  code  simulations1'  have  been  carried  out  to  investigate  the 
physical  mechanism  responsible  for  the  wave  amplification  in  the  output  section  of  a 
harmonic-multiplying  gyrotron  traveling-wave  tube.  Simulation  results  demonstrate  that 
the  injected  signal  frequency  plays  an  important  role  in  determining  the  wave  growth 
process.  This  is  because  the  initial  frequency  mismatch  as  well  as  the  resulting  electron 
energy  change  in  the  input  section  influences  the  extent  of  bunching  at  the  entrance  of  the 
output  section  and  the  slope  of  the  electron  phase  trajectory.  For  a  small  frequency 
mismatch,  the  amplification  arises  entirely  from  the  second  harmonic  generation  due  to 
the  bunching  at  the  fundamental  frequency.  As  the  mismatch  is  increased,  the  bunched 
electrons  start  to  disperse  early  in  the  output  section;  and  the  wave  amplification  is  then 
dominated  by  the  cyclotron  maser  instability.  The  electron  velocity  spread  is  observed  to 
substantially  degrade  the  device  performance. 

By  using  the  combination  of  a  waveguide  input  section  (Lj,  TE02)  that  amplifies 
the  input  signal  at  the  fundamental  cyclotron  frequency,  a  drift  section  (Lj),  and  an  output 
waveguide  section  (L0,  TE03)  that  generates  the  output  wave  at  the  second  harmonic 
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cyclotron  frequency,  the  bandwidth  of  the  phigtron18  can  be  substantially  improved.  The 
above  harmonic-multiplying  scheme  made  the  source  for  the  input  signal  more  readily 
available,  reduces  the  magnetic  weight  requirement,  and  alleviates  the  stability  problem. 

The  frequency-multiplying  amplification  process  in  the  TWA  output  section  is 
rather  involved  and  remains  to  be  addressed.  There  are  at  least  two  physical  mechanisms 
that  could  give  rise  to  the  wave  amplification  process.  The  conventional  thinking  is  that 
the  amplification  process  is  entirely  due  to  the  second  harmonic  cyclotron  maser 
instability  and  the  instability  is  believed  to  be  initiated  by  the  beam  second  harmonic 
current  component  in  the  bunched  electrons.  On  the  other  hand,  in  the  harmonic- 
multiplying  gyrotron  TWT,  the  electron  beam  is  pre-bunched  at  the  fundamental 
frequency  before  entering  the  output  waveguide.  This  fundamental  bunching  will 
continue  into  the  output  section  until  the  electrons  start  to  cross  one  another  or  the  self- 
consistent  field  grows  to  large  amplitude  to  alter  the  electron  motion.  In  a  conventional 
klystron,  the  bunched  electrons  induce  the  electromagnetic  field  in  the  output  cavity  with 
the  proper  phase  to  slow  the  beam  down  and  convert  the  beam  kinetic  energy  into 
coherent  radiation  without  invoking  any  instability.  Similarly,  the  second  harmonic  RF 
current  associated  with  the  fundamental  bunching  can  induce  an  electromagnetic  wave 
with  the  proper  phase  supported  by  the  output  waveguide  to  slow  down  the  beam.  The 
latter  mechanism  can  be  isolated  in  simulations  by  setting  the  self-consistent  field  to  zero 
and  thus  breaking  the  closed  feedback  loop  for  the  instability. 

The  parameters  used  in  simulations  for  a  35  GHz  harmonic-multiplying  gyro- 
TWT  are  Vb  =  50kV,  Ib  =  30A,  a  =  1.5,  L;  =  10cm,  Ld  =  14cm,  L0  =  20cm,  and  B0  = 
6.068kG.  In  order  to  determine  the  effect  of  varying  the  frequency  mismatch  (Aco) 
between  the  Doppler-shifted  wave  frequency  and  the  electron  cyclotron  frequency  on  the 
gain  mechanism  in  the  output  waveguide,  three  cases  with  different  input  frequency  are 
simulated.  The  spatial  evolution  of  the  second  harmonic  wave  and  the  electron  phase 
trajectories  of  30  test  electrons  in  the  output  waveguide  are  displayed  in  Fig.  5.  The 
electron  phase  trajectory  with  respects  to  the  fundamental  cyclotron  frequency  is 
evaluated  without  the  effect  of  the  RF  field  in  the  output  waveguide  according  to  <t>i(z)  = 

d>,f0)  +  T  ~r  -  IT,  where  the  argument  0  is  at  the  electron  injected  plane  and  T  is 
Y  V  ’  L  Yi(z)  Yi(0) 

the  time  elapsed.  Using  coj  =  1.03coc  (fc  is  the  cutoff  frequency  of  the  input  waveguide), 
the  initial  frequency  mismatch  Aco  is  0.01  ®c.  With  this  input  frequency,  the  wave 
amplification  (Fig.  5a)  that  takes  the  self-consistent  field  into  account  follows  very 
closely  the  result  from  the  fundamental  bunching  alone  (turn  off  the  effect  the  self- 
consistent  field  on  electron  motion).  The  electrons  remain  bunched  throughout  the  output 
waveguide  and  the  optimized  efficiency  is  about  16%.  Increasing  ®j  to  1.06coc  (Aco  = 
0.015coc),  the  optimized  efficiency  is  increased  to  22%  (Fig.  5b).  The  result  shows  that 
the  bunch  starts  to  disperse  at  about  zs  =  10  cm  and  at  about  the  same  time  the  induced 
field  from  the  fundamental  bunching  saturates.  After  zs  the  cyclotron  maser  instability 
dominates  the  gain  process.  Increasing  coj  further  to  1.08  coc  (Aco  =  0.0195coc).  The 
optimized  efficiency  reduces  to  only  8.6%  (Fig.lc).  In  this  case,  the  gain  process  is 
dominated  by  the  cyclotron  maser  instability  with  reduced  growth  rate  due  to  the  large 
beam  energy  spread  caused  by  the  dispersion  of  fundamental  bunching.  The 


14 


understanding  of  the  gain  mechanism  in  the  output  waveguide  will  have  important 
implication  for  the  design  of  future  experiments. 


z(cm) 


Figure  5.  Spatial  evolution  of  the  second  harmonic  wave  and  electron 
phase  trajectories  of  30  test  electrons  in  the  output  waveguide, 
(a)  coj  =1.03  (b)  coj  =  1.06  (c)  to,  =1.08. 
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